On quasi-exactly solvable matrix models 



R.Z. Zhdanov 

Institute of Mathematics, 
3 Tereshchenkivska Street, 252004 Kiev, Ukraine* 



Abstract 

An efficient procedure for constructing quasi-exactly solvable matrix models is 
suggested. It is based on the fact that the representation spaces of representations 
of the algebra ,3/(2, R) within the class of first-order matrix differential operators 
contain finite dimensional invariant subspaces. 

The Lie-algebraic approach to constructing quasi-exactly solvable one- dimensional sta- 
tionary Schrodinger equations as suggested by Shifman and Turbiner [jl|, ^ is based on 
the properties of the representations of the algebra s/(2, R) 

[Qo, Q±] = ±Q±, Q+] = 2Qo (1) 

by first-order differential operators. Namely, the approach in question utilizes the fact 
that the representation space of the algebra s/(2, R) having the basis elements 

Q- = -7-, Qo = x-^-n, (5+ = - 2?7.x, (2) 
ax ax ax 

where n is an arbitrary natural number, has an [n + l)-dimensional invariant subspace. 
Its basis is formed by the polynomials in x of the order not higher than n (for further 
details see the monograph by Ushveridze and references therein). 

The aim of the present paper is to extend the above Lie-algebraic approach in order 
to make it applicable to analyzing eigenvalue problems for matrix differential operators. 
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The key idea is that the basis elements of the algebra sl{2, R) are searched for within the 
class of matrix differential operators 



Q = ax)^ + V{x), (3) 
ax 

where ^{x),'r]{x) are some matrix- valued functions of the corresponding dimension. Fur- 
thermore, the representation space of sl{2, R) must contain a finite-dimensional subspace. 
Provided these requirements are met, a quasi-exactly solvable matrix model is obtained 
by composing a linear combination of the basis elements of the algebra s/(2,R) with 
constant matrix coefficients. 

Thus, to get a quasi-exactly solvable matrix model we need to solve two intermediate 
problems 

1) to solve the relations (|I]) within the class @, 

2) to pick out from the set of thus obtained realizations of the algebra si (2, C) those 
ones whose representation space contains a finite-dimensional invariant subspace. 

In a sequel we will restrict our considerations to the case when ^{x) is a scalar mul- 
tiple of the unit matrix. Given this restriction, a simple computation yields that any 
representation of s/(2,R) within the class of operators (0) is equivalent to the following 
one: 

Q- = -r^ Qo = x— + A, Q+ = x'^— + 2xA + B, (4) 
ax ax ax 

where A, B are constant N x N matrices satisfying the relation 

[A, B] = B. (5) 

Next, we have to investigate under which circumstances the representation space of 
the algebra has a finite-dimensional invariant subspace T„ with basis elements 

N 

M^) = ^ = 1, . . . , n. 

Here Ci, . . . , e/v is the orthonormal basis of the space R^. It occurs that the components 
of the functions are necessarily polynomials in x of the order not higher than n — 1. 
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Theorem 1 Let the functions fi{x), . . . , fn{x) form the basis of the invariant subspace 

^ da: ' dx 



X„ of the representation space of the Lie algebra {-i-,x4- + A), where A is a constant 



N X N matrix. Then 



0' ^ = l,...,n. (6) 

We will give a sketch of the proof for the case, when N = 2. The requirement that a 
representation space of a Lie algebra under study contains a finite-dimensional invariant 
subspace means that the following relations hold 

^/. = EA.,/„ (7) 

(x^ + a) /7 = £ L,,/; (8) 

where Aij, Lij are arbitrary complex constants, i, j = 1, . . . ,n. 
Solving (|^ yields 

N n 

/^(^) = EE (e^i,,^/^.^";' (9) 
j=i k=i 

where Cij, i = 1, . . . ,n, j = 1, . . . ,N are arbitrary complex constants and the symbol 
{A)ij stands for the {i,j)th entry of the matrix A. 
Next, from the requirement (|^) we get 

xAC + CA = e-^^Le^^a (10) 

Making use of the Cambell-HausdorfF formula and equating coefficients of the powers of 
X give the following infinite set of algebraic equations for unknown matrices L, A, C, A: 

LC = CA, (11) 
[L, A]C = AC, (12) 
{L, AyC = 0, t>2, (13) 

where 

{L, A}0 = L, {L, A}' = [{L, Ay-\ A], z > 1. 

Choosing the basis vectors ei, 6*2 in an appropriate way, we can transform the constant 
matrix A to the Jordan form. There are two possibilities 
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^ = . (15) 



Here A, Ai, A2 are arbitrary constants. 

Case 1. Let the matrix A be of the form (0). If we denote the first and the second 
columns of the 2x2 matrix C as Ci and C2, then equations (^) can be rewritten to 
become 

LCi = ACi, LC2 ~ Ci + AC2- 
Next, using equations (|12[) , (|13D we obtain the following relations 



LA^Ci = XjA^Ci, (16) 
LA^Ca = A^Ci + XjA^C2- (17) 



Here Aq = A, Aj+i = Aj + 1, j = 0, 1, . . .. 



It follows from (|l^) that aj+i = A'Ci, i > are eigenvectors of the n x n matrix 
L corresponding to eigenvalues Aj = A + i and, what is more, since these eigenvectors 
correspond to distinct eigenvalues, they are linearly independent. As there are at most n 
linearly independent eigenvectors of the matrix L, the relation 

A'^Ci = (18) 

with some m < n holds true. 

Combining (|T^ and ([T^) yields 

LA"^+^Ci = A^+,A"^+^C2, I > 0. (19) 

Hence, we conclude that the vectors Sm+i = A'"+*C2, i > are eigenvectors of the matrix 
L forming together with the vectors the system of its linearly independent 

eigenvectors. As an n x n matrix has at most n linearly independent eigenvectors, the 
relation 

A'^(72 = (20) 

holds true. 

In view of (|T^), ( PUj ) the matrix C = (Ci, C2) satisfy the following matrix equation 

A"C7 = 0. 
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Due to this fact, reads as 

N n / 2 n-l \ 

M^) = E E 1 + + -A^ + ■ ■ ■ + T-^rryA"-^ C,,e,. (21) 
j=i k=i \ ^- / ik 

Thus, the components of the vectors fj are polynomials of the order not higher than n — 1, 
which is the same as what was to be proved. 

Case 2. We turn now to the case when the matrix A is given by (|l^). Denoting the first 
and the second columns of the 2x2 matrix C as Ci and C2, we rewrite equations (0) as 
follows 

LCi = AiCi, LC2 = A2C2- 
Next, using equations (|T^, (|T^ we obtain the relations: 

LA^Ci = ajA^Ci, (22) 
LA^C2 = l3jA^C2, (23) 

where 

tto = Ai, Oj+i = ttj + 1, j = 0,l,..., 
/?o = A2, Pj+i = (3j + l, J = 0,1,.... 

Thus, the vectors Cj+i = A*Ci, i = 0, 1, . . . are eigenvectors of the n x n matrix L 
with eigenvalues = Ai + i, i = 0, 1, . . .. As these eigenvalues are distinct, the vectors 
Cj are linearly independent. Taking into account that an n x n matrix can have at most 
n linearly independent eigenvectors we conclude that A"Ci = 0. Similarly, we get the 
relation A"C2 = 0. Hence it follows that the expression (|^) takes the form (^Tj), which is 
the same as what was to be proved. • 

Consequently, the most general finite-dimensional invariant subspace of the represen- 
tation space of the algebra (|^) is spanned by vectors whose components are finite-order 
polynomials in x. We postpone a detailed study of these representations with arbitrary 
N for future publications and concentrate on the case N = 2. There are two families 
of inequivalent finite dimensional representations of the algebra s/(2, R) having the basis 
elements (^ 
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Here n, m are abritrary natural numbers with n> m. 

Representations of the form (P, (|2^) are the direct sums of two irreducible represen- 
tations realized on the representation spaces 

7^l = (ei,xei, . . . ,x"ei), 7^2 = {e2,xe2, . . • ,x"'e2), 

where ei = (1,0)^, 63 = (0, 1)^. 

Next, representations (|), (^5]) are also the direct sums of two irreducible representa- 
tions realized on the representation spaces 

TZi = {nil, • ■ • 5 nx^ei + jx^~^e2, . . . , nx^ei + nx^~^e2), 
K2 = (e2,xe2, . . . ,x""^e2). 

We will finish the paper with an example of utilizing the above results for obtaining an 
exactly solvable two-component Dirac-type equation, which is one of the two differential 
equations composing the Lax pair for the cubic Schrodinger equation (see, e.g. @, H). 
Consider the following two-component matrix model: 

Hw = i(axa2 + bai)—— + (ciai + (c2 + ia) (12) w = Xw, (26) 
ax 

where a,b,Ci,C2 are arbitrary real parameters with ab and cri,o"2 ai'e 2x2 Pauli 
matrices. As 

H = ibaiQ^ + iaa2Qo + c\0\ 020-2, 

where Q-,Q+ are given by (^) with A = 1, the operator Ti transforms the {2n + 2)- 
dimensional vector space 

fj{x) = x^'^ei, fn+j+iix) = x^~^e2, j = l,...,n + l 

into itself. This means that there exists the constant {2n -|- 2) x (2n -|- 2) matrix H such 
that 

2n+2 

n.n= ^ H,,fl i = l,...,2n + 2. 



Hence it immediately follows that the vector-function 

2n+2 

E 
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is the solution of the system of ordinary differential equations (|26| ) with A = Aq, provided 
(«!, . . . , a2n+2) is an eigenvector of the matrix H with the eigenvalue Aq. 
Making a transformation 

X = —sinh^ay), 
a 



wix) 



cosh(a?/))^^^ exp < — (ci arctan sinh(ay) + C2 In cosh(a?/)) 



z 
a 

X exp{— icTs arctan sinh (ay)} ^^(2/) 



we reduce (p6|) to the Dirac-type equation 

^^i^ + cr2n2/)^(l/) = AV^, (27) 

where 

N a^C2 — fo^Ci sinh(ay) 

= — T:r~\ — 

aocosh[ay) 

is the well-known hyperbolic Poschel- Teller potential. It has exact solutions of the form 

ip{y) = (cosh(a?/))^^^^ exp I - (ci arctansinh(a?/) + C2 In cosh(ay)) 

I a 

2n+2 ^ \ 

X exp {z arctan sinh( ay) as} ^ ajfj -sinh(az/) . 

j=i \" / 

As the potential V does not depend explicitly on n, the order of the polynomials P„, Qn 
may be arbitrarily large. This means that the Dirac equation (^) with the hyperbolic 
Poschel- Teller potential is exactly-solvable. 

In we suggest an alternative approach to construction of quasi-exactly solvable 
stationary Schrodinger equations based on their conditional symmetry. We believe that 
a similar idea should work for matrix models as well. It is intended to devote one of 
the future publications to a comparison of the conditional symmetry and Lie-algebraic 
approaches to constructing quasi-exactly solvable Dirac-type equations (p^). 
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